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We  consider  the  theory  of  chemisorption  on  stepped  surfaces  for  systems 
for  which  a  lattice  gas  in  thermodynamic  equilibrium  is  an  appropriate  model. 
Using  statistical  mechanics  we  demonstrate  that  by  comparing  LEED  (Low  Energy 
Electron  Diffraction)  results  for  flat  and  stepped  substrates  one  can  determine 
the  change  in  adsorption  energy  at  terrace  edges.  Thus  one  can  answer  the 
important  question  of  whether  adatoms  bind  more  or  less  strongly  at  this  kind 
of  defect  site.  We  also  examine  some  new  and  interesting  features  of  LEED 
spot  shapes  and  integrated  intensity  behavior  for  adsorption  on  stepped  sub¬ 
strates  when  the  overlayer  is  partially  disordered.  In  a  companion  paper,  we 
consider  0  on  a  certain  stepped  W(110)  surface.  We  demonstrate  there,  via  Monte 
Carlo  calculations,  that  the  binding  energy  is  less  strong  on  either  terrace 
edge  for  this  particular  adsorbate,  surface  and  defect. 

I.  Introduction 

The  detailed  behavior  of  adsorption  at  steps  is  a  question  of  fundamental 
importance  to  surface  physics.  The  growth  kinetics  and  equilibrium  shape  of 
crystals  can  be  determined  by  whether  adsorption  occurs  at  step  edges  or  on 
the  flat  surface.  Overlayer  ordering  can  be  strongly  influenced  by  terraces 
if  domains  of  different  azimuthal  ordering  are  crystallographically  equivalent 
on  the  flat  surface  but  distinguished  by  step  edge  direction.  Reactions  in 
heterogeneous  catalysis  are  often  associated  with  "active  sites"  at  defects,  and 
the  reaction  kinetics  may  be  profoundly  influenced  by  adsorption  energies  at 
those  sites. 
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In  this  paper  we  present  a  detailed  theory  of  LEED  results  for  certain 
chemisorbed  systems  on  stepped  surfaces.  We  also  show  how  information  about 
changes  in  adsorption  energies  at  terrace  edge  sites  may  be  extracted  from 
those  results. 

Section  II  presents  a  general  theory  of  LEED  scattering  from  overlayers 
on  stepped  surfaces,  including  the  effects  of  disorder  in  the  overlayer  when  a 
lattice  gas  model  is  appropriate.  We  consider  the  overlayer  spot  splitting  and 
integrated  intensity  in  the  kinematic  (single-scattering)  approximation  and  also 
point  out  that  qualitatively  new  LEED  features  may  appear  at  low  or  high  coverages. 

Our  general  conclusions  are  illustrated  with  a  model  adsorption  system 
calculation  at  T  =  0  in  Section  III.  We  consider  the  effects  of  the  change  in 
edge  site  adsorption  energy  on  overlayer  LEED  spot  splitting  and  intensity  as 
a  function  of  coverage. 

In  a  companion  paper1,  we  apply  the  theory  developed  here  to  0  on  a  certain 
stepped  W(110)  surface.  Using  experimental  results  and  existing  models  for  the 
stepped  and  flat  surface  adsorption  systems  as  input  in  a  Monte  Carlo  calculation, 
we  show  that  the  edge  site  adsorption  energy  is  less  strong  at  either  terrace 
edge  than  for  the  flat  surface. 

The  lattice  gas  theory  employed  here  makes  use  of  the  properties  of  thermo¬ 
dynamic  equilibrium.  This  has  some  important  advantages  in  surface  systems  that 
have  not  been  fully  exploited,  especially  in  chemisorption  studies.  First, 
properly  chosen  systems  are  characterized  by  only  a  few  parameters.  Since 
there  exist  several  well -developed  methods  for  calculating  thermodynamic 
quantities,  one  can  determine  these  parameters  with  considerable  confidence. 

For  instance,  in  an  adsorption  system  many  observable  effects,  such  as  LEED 
results,  depend  only  on  the  adsorption  energy  and  adatom-adatom  (AA)  Interactions. 
This  Is  a  consequence  of  the  general  fact  that  all  kinetic  effects  are  Irrelevant 
to  equilibrium  properties.  This  situation  Is  especially  advantageous  when  one 
Is  dealing  with  a  lattice  gas  model  with  finite  range  AA  forces.  There  are  only 
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a  few  energy  parameters  in  the  problem,  since  the  adatoms  sit  at  discrete 
sites  and  are  only  influenced  by  a  few  near  neighbors.  Second,  although 
equilibrium  studies  certainly  cannot  solve  all  the  important  questions  in 
surface  physics,  they  can  provide  a  secure  starting  place  from  which  to  attack 
other  more  difficult  problems  in  complicated  systems.  One  such  application  of 
equilibrium  studies  on  a  flat  surface  is  demonstrated  in  the  companion  paper 
on  0  on  stepped  W(110),  and  we  believe  it  will  not  be  the  last. 

II.  LEED  Scattering  from  Overlayers  on  Stepped  Surfaces 

In  this  section  we  present  a  theory  of  LEED  scattering  from  overlayers 
on  stepped  surfaces.  There  have  been  several  treatments  of  LEED  scattering 
from  stepped  surfaces  ,  but  to  our  knowledge  none  of  them  have  considered  the 
effects  due  to  partial  ordering  of  the  adlayer  that  we  present  below.  Our 
treatment  uses  the  single  scattering  (kinematic)  approximation  throughout. 

While  multiple-scattering  (dynamic)  effects  constitute  a  very  important  part 
of  the  LEED  signal,  in  certain  cases  they  may  be  ignored  when  considering  an 
overlayer  beam  and  a  small  angular  region  around  it.  Previous  work  has  shown 
that  overlayer  spot  shapes  and  trends  in  the  intensity,  rather  than  their 
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absolute  numerical  values  ,  are  well  approximated  by  the  single  scattering 
approximation. 

Now  consider  a  stepped  surface.  We  assume  that  this  surface  is  composed 
of  identical  terraces,  with  equivalent  points  on  neighboring  terraces  separated 
by  a  fixed  vector  a,  and  that  there  are  N$  such  terraces  within  the  coherence 
length  C  of  the  LEED  apparatus.  If  Jc  Is  the  momentum  transfer  (scattering) 
vector  Xt  a  kf  -  JSq*  =  incoming  (scattered)  electron  wave  vector, 

IkJ  *  (Jc^l )  then  the  total  scattering  amplitude  A(k)  is 
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Here  A.(Jc)  is  the  amplitude  due  to  the  jth  terrace,  and  £  is  a  (real  space) 

J 

vector  connecting  identical  points  on  adjacent  terraces.  For  monatomic  height 
steps  and  uniform  terraces  £  ~  Ngg.  If  the  terraces  are  identical  and  we  use 
the  single  scattering  approximation 

j(t)  t  (c  4 


In  Eq  (2)  we  have  made  the  important  assumption  that  the  overlayer  can  be 
described  by  a  lattice  gas  model.  This  means  the  adatoms  on  a  given  terrace 
sit  on  a  fixed  set  of  N  equivalent  sites  &,  with  occupation  variable  n^  =  0(1) 
for  an  unoccupied  (occupied)  site.  These  sites  are  at  positions  with  respect 
to  the  origin  of  £.  (Mote  that  the  location  of  these  adatom  sites  with  respect 
to  the  substrate  is  irrelevant  in  the  lattice  gas  picture).  Some  of  the  other 

4 

assumptions  involved  in  using  a  lattice  gas  model  are  reviewed  by  Domany  et  al  . 

In  writing  Eq.  (2)  we  have  ignored  an  overall  multiplicative  factor  -  the  adatom 
scattering  cross-section  -  which  is  usually  sufficiently  slowly  varying  with 
scattering  angle  to  be  taken  constant,  in  the  vicinity  of  any  diffraction  feature. 
Now  the  LEED  intensity  for  the  beam  k^  is  just 


N* 


(3) 


XU)*  <IAU)F>  -  H 

d.*  - 1 

In  Eq.  (3)  the  angular  brackets  denote  an  average  over  adlayer  configurations. 

We  assume  from  this  point  on  that  the  adlayer  is  in  thermodynamic  equilibrium, 
which  means  that  at  a  given  temperature  T  and  coverage  0  the  average  will  be 
determined  by  adatom-adatom  (AA)  forces  and  any  variation  In  the  adatom  adsorption 
energy  e  (AS  force)  due  to  the  steps.  Now  in  the  case  of  strong  chemisorption, 
e  Is  generally  several  eV  (for  0/W,  e  -  5eV).  Further,  e  is  mainly  determined 
by  the  few  substrate  atoms  surrounding  an  adsorption  site.  Since  the  configuration 
of  neighbors  Is  different  for  a  site  just  below  or  just  above  a  step,  one  expects 
a  change  6c  In  adsorption  energy  there  on  the  order  of  one  eV,  so  that  6e  »  kT. 


Hence  terrace  edge  sites  for  which  6c  >  0  will  be  the  last  to  be  occupied  as  the 
coverage  0  is  increased.  On  the  other  hand,  sites  for  which  6c  <  0  will  be 
occupied  first,  and  can  serve  as  nucleation  sites  for  overlayer  domains.  This 
distinction  can  have  a  strong  effect  on  the  0  dependence  of  the  LEED  scattering, 
which  we  examine  in  detail  below.  Such  effects  can  then  be  used  to  determine 
the  sign  of  6e  from  LEED  results  for  a  given  stepped  chemisorption  system. 

Now  for  many  important  adsorption  systems  the  AA  interactions  are  short 
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ranged  .  Since  the  adsorption  sites  immediately  on  either  side  of  a  step  will 
be  either  occupied  or  empty  according  to  whether  6e  <  0,  the  only  AA  interactions 
that  can  "cross"  a  step  and  influence  the  next  terrace  will  be  long  range  and 
hence  quite  weak.  Thus  we  assume  that  the  statistical  average  in  Eq(3)  can  be 
taken  independently  on  each  terrace.  Hence 

I(fe}  r  KU  <lf(U|V  *£  (</a>>r«pUk-i  <]'-*>)  C</> 

r  Ns(  -  I  <X^)>r  )  v  InMj  fe) 

In  Eqs.  (4)  and  (5)  the  statistical  average  refers  to  a  single  terrace,  and 
the  interterrace  interference  Is  expressed  via 

h(Ns,*)=£S  **(>(<: -  Sin  C^X/r)  /sinft/l)  &  ) 
This  function  is  familiar  from  the  simple  diffraction  grating.  For  N$  >  1  It 

O 

has  peaks  of  height  N  and  width  27r/N  centered  at  x  *  0,  ±  2ir,  ±  4wt  and  Is 
quite  small  in  between.  The  peak  area  rapidly  approaches  2ttNs  as  N$  grows. 

Eq  (5)  Is  a  central  result  of  our  theory.  Its  first  term  refers  to  a 
single  terrace  and  is  non- zero  only  if  there  is  some  degree  of  statistical 
disorder  of  the  adatoms.  This  term  contains  some  new  physics  that  to  our 
knowledge  has  not  been  examined  previously.  The  second  term  In  Eq  (5)  Includes 
the  effects  of  Interterrace  Interference  via  the  function  h.  If  one  has  a 
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2  2 

"complete"  (perfect)  adlayer  then  <|p|  >  =  j<p>|  and  the  peaks  in  h  express  the 

condition  of  constructive  Interference  between  steps.  This  gives  rise  to  a 

characteristic  splitting  of  the  LEED  spots  in  the  £  direction  from  which  the 

2 

average  step  width  and  height  may  be  determined  .  If  there  is  disorder  in 
the  adlayer,  the  situation  is  more  complicated.  In  general,  the  second  or 
"split"  term  in  Eq(5)  will  be  non-zero.  In  addition  there  will  also  be  a  con¬ 
tribution  to  the  LEED  intensity  from  the  first  term,  due  to  the  disorder,  that 
refers  to  one  terrace  only  and  is  not  split.  In  what  follows  we  consider  the 
influence  of  the  two  terms  in  Eq(5)  on  overlayer  LEED  spot  splitting  and  integrated 
intensity.  As  we  will  see,  the  sign  of  6e  (the  change  in  adsorption  energy  at 
terrace  edge  sites)  plays  an  Important  role  in  these  effects. 

It  is  interesting  to  compare  Eq(5)  with  the  van  Hove  expression  for  thermal 
neutron  scattering.6  The  first  term  in  Eq(5)  corresponds  to  the  incoherent 
neutron  scattering  cross-section  and  p  corresponds  to  the  (spin  dependent)  neutron- 
single  nucleus  scattering  length. 

We  now  consider  uniform  terraces  that  are  long  compared  to  their  width. 

Then  we  need  only  consider  two  Se,  on  either  long  edge.  In  general,  for  small 

widths  and  attractive  6e  on  one  or  both  edges,  one  would  expect  most  adatoms  to 

be  "frozen"  In  a  single  ordered  configuration  as  long  as  0  and  T  are  not  too 

large  (l.e.  for  T  <  T  ,  where  T  *  ordering  temperature  on  the  flat  surface  and 
**  c  c 

9  <  0C  where  Qq  is  the  coverage  of  a  "complete"  ordered  overlayer  on  the  flat 

surface).  In  this  case  the  second  term  in  Eq(5)  will  be  more  important 

and  splitting  may  be  observed.  The  first  term,  which  is  a  measure  of  adlayer 

disorder,  will  in  general  be  relatively  small.  If  both  edges  are  repulsive, 

the  situation  is  more  complicated.  Then  for  9  <  9C,  one  can  have  "floating" 

2  - 

regions  of  ordered  adlayer,  with  a  randomn  phase  relation  so  that  |<pfc>|  *  0 
and  there  will  be  little  or  no  splitting.  As  0  grows,  there  are  two  possibilities. 

If  the  terrace  width  Is  commensurate  with  an  Integral  number  of  adlayer  unit 

2  2 

meshes  the  overlayer  "locks  in"  to  a  single  ordered  arrangement  so  that  |<p>|  *  <|p|  > 


and  the  spot  is  split.  If  the  width  is  incommensurate  with  an  integral  number 

a 

of  adlayer  unit  meshes,  one  can  have  { <p> |  =  0  even  for  6  =  9C.  On  a  real 
surface  one  expects  a  distribution  of  terrace  widths  so  that  a  combination  of 
both  these  possibilities  will  occur.  These  various  cases  are  illustrated  by 
our  model  in  Section  III. 

It  is  also  instructive  to  consider  the  limit  of  Eq ( 5)  for  infinite  terrace 

width.  In  this  case  the  number  of  sites  N  »,  and  the  number  of  terraces  M$  1 

as  the  width  approaches  the  LEED  coherence  length  £.  When  this  happens,  the 

interterrace  interference  function  h  -*•  1  and  the  right  hand  side  of  Eq(5)  reduces 
2 

to  <  | p |  >,  which  is  the  proper  limit.  For  N  -*•  »,  this  function  will  be 

independent  of  6e  except  in  certain  special  cases.  Note  also  that  as  N  grows, 

2  9  2  9  2  2 
<  | p |  >  -*■  0  (N  ),|<p>|  -►  0(N  ),  and<  |p|  >  -  (<p>|  a  NTx,  where  x  is  the 

subceptibility  of  <p(k)>  to  a  field  coupling  to  p(k).  One  has  x  *  0(1),  but 

X  may  become  large  near  an  incipient  phase  transition  (a  true  phase  transition 

cannot  occur  in  a  two  dimensional  system  with  finite  range  forces).  Hence  for  N 

large  but  not  infinite  the  first  term  in  Eq(5)  may  become  important. 

Because  of  finite  instrumental  response  function  effects,  an  important 

quantity  for  comparison  with  LEED  experiments  is  the  integrated  intensity 


4  =  $$*(*)  4M> 


where  a  and  b  are  unit  mesh  vectors  and  the  integral  extends  over  an  area  a  in 
reciprocal  space  around  the  overlayer  diffraction  beam,  a  is  defined  by  in¬ 
strumental  parameters  and  beam  energy.  It  is  apparent  from  Eq(5)  that  even  if 
2 

| <p> |  is  sizable  Its  contribution  to  I  may  be  relatively  small,  depending  on  a 

2 

and  the  shape  of  |<p(k)>|  as  a  function  of  k.  For  Instance,  suppose 
2 

| <p> |  s  constant.  Then  since 


ii 


there  is  a  tendency  for  the  two  contributions  of  |<p>|  to  *'  to  cancel.  To 

a 

understand  this  point  in  more  detail  let  define  an  adlayer  beam  and  i  be 

2 

a  unit  vector  parallel  to  the  projection  of  £  onto  a  terrace.  Then,  if  j<p(Jk)>|  = 

A 

constant  for  jc  in  the  range  +  ti  i  and  if  there  is  a  maximum  of  h  at  x  =  the 

2  >  9 

contribution  of  | <p> |  to  will  be  small  if  the  integrated  intensity  is 

measured  over  this  range  of  J<  values.  Notice  that  for  N$  reasonably  large 

(>2  or  3)  and  for  sufficiently  wide  terraces,  the  peaks  in  h  will  be  quite 

narrow  (Ak  -  ).  Hence  h  will  act  like  a  delta  function  in  the  £  direction 

and  its  affect  on  ~  will  be  to  single  out  one  k^  value.  Further,  one  expects 
2 

|<p(k)>|  to  "spread  out"  as  a  function  of  k^  for  g  small  enough.  Hence  it 

may  happen  that  the  LEED  instrument  response  function  is  such  that  the  net 

2 

contribution  of  |<p>|  to  is  small.  This  in  fact  occurs  in  the  model  calcula¬ 
tion  for  0/stepped  W(110)  reported  in  the  companion  paper J 

There  is  an  interesting  possibility  of  new  LEED  features  appearing  in  the 
stepped  chemisorption  systems  we  are  considering.  Suppose  Se  <  o  for  a  certain 
edge.  Then  the  sites  along  that  edge  will  fill  up  first,  while  all  others  on 
the  terrace  remain  empty.  Since  the  substrate  configuration  is  quite  different 
at  the  edge  than  elsewhere  in  the  terrace,  one  would  expect  the  adatom-adatom 
(AA)  interactions  to  be  different  as  well.  If  the  AA  edge  energies  favor  a 
particular  one-dimensional  order  and  the  temperature  is  not  too  large,  new  LEED 
spots  will  appear.  If  one  edge  has  6e  <  o  and  the  terrace  is  D  sites  wide, 
these  spots  will  reach  a  maximum  intensity  for  some  coverage  o  <  9  <  ^  .  The  spots 
will  be  streaked  in  the  direction  perpendicular  to  the  step  edge,  since  they 
are  due  to  one-dimensional  order.  There  should  also  be  a  spot  splitting  in 
this  direction  due  to  the  fixed  phase  relation  between  edge  sites.  If  the 
terraces  are  randomnly  distributed  In  the  direction  parallel  to  the  step  edge, 
there  will  be  no  splitting  In  that  direction.  For  a  repulsive  edge  (<5e>o)  similar 
effects  can  occur  for  coverages  1  -  ^  <  e  <  1 .  Note  that  the  two  terrace  edges 
are  different,  hence  one  can  have  two  new  sets  of  AA  interactions.  The 


corresponding  LEED  features  will  occur  simultaneously  at  low  9  if  both  <5e<o, 
at  high  9  if  both  6e>o,  or  separately  at  low  and  high  9,  if  6e<o  for  one  edge 
and  6e>o  for  the  other.  To  our  knowledge  such  LEED  features  have  neither  been 
predicted  nor  observed  previously,  and  we  urge  that  experiments  be  done  to 
confirm  these  ideas. 

III.  Model  Calculation 

We  illustrate  the  general  points  made  above  with  a  model  calculation  at 
T  =  0.  Consider  a  square  lattice  of  adsorption  sites  with  terraces  of  width 
D  sites  in  the  [1,0]  direction  and  length  L  sites  in  the  [0,1]  direction. 

Suppose  the  adatom-adatom  (AA)  interactions  are  such  that  (2X1)  adlayer  ordering 
occurs.  The  "perfect"  (complete)  overlayer  has  an  adatom  occupying  every 
other  site  in  the  [1,0]  direction,  and  every  filled  site  is  part  of  a  filled 
row  in  the  [0,1]  direction.  A  simple  (but  not  unique))  set  of  AA  interactions 
giving  rise  to  this  type  of  order  is  nearest  neighbor  attraction  on  the  [0,1] 
direction  along  with  nearest  neighbor  repulsion  and  second  neighbor  attraction 
on  the  [1,0]  direction. 

We  restrict  ourselves  to  T  =  0  to  be  able  to  calculate  analytically.  In 
general,  we  expect  the  conclusions  we  draw  at  T  =  o  will  be  valid  for  a  con¬ 
siderable  range  of  finite  temperatures.  This  expectation  is  born  out  by  the 
results  for  our  finite  temperature  calculation  for  0  on  stepped  W(110)  in  the 
companion  paper\  There  may  well  be  more  subtle  effects  near  incipient  phase 
transitions,  but  we  are  not  concerned  with  them  here. 

We  now  consider  several  aspects  of  the  LEED  scattering  predicted  by  our 

model  at  T  =  0.  The  features  of  interest  are  (i)  the  behavior  of  the  quantities 
2  2 

< | p |  >  and  | <p> 1  and  the  LEED  Intensity  I  vs.  coverage  0  for  the  overlayer  beam, 

(ii)  LEED  spot  splitting  effects  and  (111)  the  relative  size  of  the  contributions 

2  2  1 
of  < | p (  >  and  | <p>|  (cf.  Eq(5))  to  the  Integrated  LEED  intensity  at  certain  8 

values.  There  are  three  possibilities  for  the  change  in  adsorption  energy  at 


terrace  edge  sites:  (a)  both  6e>o  (both  edges  less  attractive  than  the  terrace) 
which  we  denote  ++,  (b)  both  <Sc<o  (both  edges  more  attractive  than  the  terrace, 
denoted  — ,  and  (c)  one  edge  with  6e>o  and  one  with  6e<o,  denoted  +-.  Further 
we  distinguish  in  each  case  whether  D  is  even  or  odd  since  this  difference  can 
be  significant.  We  make  no  assumptions  here  about  AA  interactions  along  the 
edges  and  therefore  do  not  consider  the  possibility  of  new  LEED  spots  mentioned 
at  the  end  of  Section  II. 

Now  let  a^  (bj  be  a  vector  between  neighboring  adsorption  sites  in  the 
[1,0]  (  [0,1]  )  direction.  A  (2x1)  overlayer  beam  may  be  defined  by  j^,  where 
^•a  =  it,  k^-b  =  0.  If  the  experiment  is  set  up  so  that  =  (2  n+l)TT,  n  =  integer, 
the  second  term  in  Eq(5)  will  not  contribute  to  the  intensity  of  the  overlayer 
beam  I (kg)'  When  considering  splitting  we  will  also  be  interested  in  a  nearby 
beam  k^  chosen  to  coincide  with  a  maximum  of  h  in  Eq(5),  kjj  =  k^  +  Ak/£  =  ±ir. 

There  is  a  third  overlayer  beam  defined  by  k_  =  1  that  is  of  interest;  for  this 

'  o 

beam  (N  .  k  ' -q)  =  N,.  One  can  always  find  such  a  k  vector  since  o<h<Nc  .  This 

S  O  5  - ^ 

beam  is  useful  in  discussing  the  intensity  I  since  by  Eq(5) 

1  (L<o)  -  M*  =■  Ms  WjCWlN  (i) 

The  last  equality  holds  because  the  statistical  average  is  independent  of  the 

component  of  J<  perpendicular  to  the  terrace. 

Mow  we  consider  case  (a),  ++,  for  D  =  odd.  At  low  0,  there  is  considerable 

2 

cancellation  in  |<p(kQ)>|  »  since  the  (2x1)  domains  may  nucleate  anywhere  except 

at  an  edge  site.  The  maximum  value  of  this  quantity  comes  at  0m  =  1/2  <  1/2. 

2  2 

< I p( i^) !  >  will  be  larger  than  J <p> |  at  low  0,  but  will  also  reach  a  maximum  with 

2 

the  same  value  as  |<p(k0)>|  at  0m.  This  may  be  seen  easily  by  finding  the  state(s) 
of  minimum  energy  at  each  0.  Note  that  the  maximum  LEED  intensity  for  a  flat 
surface  (2x1)  overlayer  comes  at  0  =  1/2.  Now  will  go  to  zero  at  0m>  but 


I(J^')  will  have  a  maximum  there.  Since  0^  =  1/2  -  jq  for  the  ++ 

D  =  odd  case,  the  maximum  in  I(JSq')  1s  shifted  to  lower  coverage  values  than 
for  the  flat  surface.  This  is  simply  a  consequence  of  the  fact  that  less 
adsorption  sites  are  available,  which  implies  a  "complete"  adlayer  forms  at 
a  lower  6  value.  By  our  arguments  in  Section  II,  one  can  expect  similar  results 

for  the  integrated  intensity  -if  the  integration  area  is  properly  chosen, 

2 

since  the  contribution  of  | <p> |  may  be  quite  small.  This  is  confirmed  in  the 
companion  paper^. 


For  0=0  one  has 
m 


<|c\S*  UpO-L.1"  stu1-  (d-iJ  k -ft)  /  If  .a-) 


0  oj 


Hence,  using  g  =  Da,  one  finds 

X(W.)  =  L"  7>Vtt"  (,() 

Since  I(J^)  =  0,  the  overlayer  spot  is  strongly  split  (for  a  real  system  various 
effects  that  we  are  ignoring  such  as  multiple  and  diffuse  scattering  will  make 
I(k  )  >  0).  If  one  calculates  the  integrated  intensity 

4  =  f  ruo  OO 

J  xTT 

-n  -  iT/j) 

2 

where  the  integral  is  along  a  line  including  and  Jc ^ ,  one  finds  that  the  <|p|  > 

term  in  Eq(5)  contributes  about  5/4  of  the  total  integrated  intensity  for  D  =•  It. 

2 

For  case  (a),  ++  and  D  =  even  the  behavior  of  <|p(k^)|  >is  similar  to 

D-2 

that  above,  except  that  its  maximum  is  at  0m  =  1/2  -p-  ,  further  below  1/2 

2 

than  for  the  D  =  odd  case.  I <p(J<0 )> I  »  however,  goes  to  zero  at  0m  due  to 
phase  cancellation  effects.  Thus  both  I(k^)  and  I(k^')  show  maximum  at  0  =  ©m. 

At  this  coverage,  one  finds 

<|jT>  =  Ll  sin  (i-OCfe-frV  /  03/ 


i<j»'  * 


o*/) 


cos'  'll  Ik'*) 


Hence 


T(k0  )  r  Nl  b  /i  )  a  ^ 

r  Us,)  =  (  l  -  Or/i))1 }  Ci/TrJ1"  4  ^  ^/i)  CVtt}L 

-  (0/jr)L-f)  +  K/Sv  O^J 

so  that 

x(k,)/xCU  ~  =  ~>/X 

0>-i  /  O 

if  we  take  D  =  10,  Ng  =  5.  Therefore  the  overlayer  spot  is  not  split.  De- 

2 

fining  •  as  in  Eq ( 1 2)  one  finds  that  the  contribution  of  the  | <p> J  term  in  I ( k) 

2 

is  about  10%  and  the  j <p> |  h  term  about  2%  of  the  total. 

For  case  (b),  —  and  D  =  even,  there  is  only  one  overlayer  arrangement 

2  2 

possible  at  low  T  and  < ( p(Jl)  (  >  =  (<p(kj>[  for  all  k.  Also,  < f P )  f  >is  very 

small  due  to  phase  cancellation  effects.  For  D  =  odd,  the  maximum  in  either 

term  comes  at  0m  *  1/2  -g—  =  1/2  +  .  On  a  real  surface  with  both  6e<o  one 

would  expect  a  mixture  of  odd  and  even  width  terraces.  Hence  the  maximum  in 

Ifk^)  or  will  occur  for  coverages  aboye  the  flat  surface  value  of  1/2. 

For  case  (c),  +-,  there  is  again  only  one  overlayer  arrangement  possible, 

2  2 

so  < | p |  >  =  | <p> |  .  For  D  =  even  the  maximum  comes  at  em  =  1/2,  for  D  =  odd 

«  e™  ■  W  T1  ■  '/2  -  Id  • 

Notice  that  one  can  distinguish  between  the  three  cases,  ++,  —  and  +- 
(both  edges  less  attractive,  more  attractive,  or  one  of  each)  by  simply  looking 
at  the  coverage  0m  at  which  the  observable  quantity  1(1^')  peaks  vs.  the  flat 
surface  value  (0  *  1/2  here).  Assuming  an  even  mixture  of  odd  and  even  terrace 
widths  one  finds 


(S 


For  D  =  10  which  is  achievable  by  cutting  crystals  at  an  angle,  the  difference 
between  these  three  cases  is  easily  within  the  resolution  of  Auger  spectroscopy. 
Thus  one  can  determine  whether  terrace  edge  sites  adsorption  energies  are 
greater  or  less  than  on  the  terrace  by  looking  at  changes  in  the  coverage  de¬ 
pendence  of  appropriately  chosen  LEED  features.  Note  that  9m  is  respectively 
less  than,  greater  than  or  close  to  the  flat  surface  value  depending  on  whether 
there  are  effectively  fewer,  more,  or  the  same  number  of  adsorption  sites 
available  on  the  stepped  surface  than  the  flat  surface.  Thus  we  expect  the 
shift  induced  by  the  changes  in  adsorption  energy  on  the  stepped  surface  to  be 
a  very  general  effect,  determined  only  by  the  sign  of  the  two  6c  and  not 
specific  to  the  model  calculation  presented  here.  Further  consequences  of  this 
argument  for  theory  and  experiment  are  explored  in  the  companion  paper^ . 
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